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Abstract 



A gauge-invariant formulation of Fermi's Golden rule is proposed. We 
shall rivisit the conventional description of carrier-phonon scattering in the 
presence of high electric fields by means of a gauge-invariant density-matrix 
approach. We show that the so-called intracollisional field effect — as usually 
accounted for — does not exist: it is simply an artifact due to the neglect of 
the time variation of the basis states which, in turn, leads to a ill-defined 
Markov limit in the carrier-phonon interaction process. This may account for 
the surprisingly good agreement between semiclassical and rigorous quantum- 
transport calculations. 
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Since the early days of quantum mechanics |IJ the field-induced coherent dynamics of an 
electron wavepacket within a crystal, known as Block oscillations (BO), has attracted signif- 
icant and increasing interest ||. Indeed, the problem of properly describing the scattering- 
free motion of an electron in a solid has led to a three-decade controversy on the existence 
of BO ||; This originated from the different approaches employed for the description of 
the applied field, namely the vector-potential or accelerated- Block- state picture {§] and the 
scalar-potential or Wannier-Stark description ||. As discussed in ||, these two pictures are 
now recognized to be fully equivalent, since they correspond to different quantum- mechanical 
representations connected by a gauge transformation. 

The presence of scattering as well as tunneling processes strongly modifies such ideal BO 
scenario j7|. In particular, non-elastic interaction mechanisms — like carrier-LO phonon 
scattering — tend to spoil such coherent dynamics, leading to a nearly semiclassical or 
Boltzmann-like transport picture. In the presence of strong electric fields, however, the 
use of the conventional scattering picture — involving transitions between field-free Bloch 
states within Fermi's golden rule — becomes questionable. 

As originally pointed out by Levinson j§| and by Barker and Ferry ||, the effect of the 
field during the scattering process, usually referred to as intracollisional field effect (ICFE), 
may lead to significant deviations from the semiclassical scenario. On the one hand, the 
role played by the ICFE has been extensively investigated by means of rigorous quantum- 
transport approaches [pT0| P~3T - Their application, however, was often limited to highly sim- 



plified physical models and conditions, thus preventing from any quantitative comparison 
with experiments. On the other hand, strong effort has been devoted to incorporate the 



ICFE within conventional — and more realistic — Monte Carlo simulations []14j] . In this case, 
the basic idea is that, due to the field-induced carrier drift, energy conservation in the scat- 
tering process is relaxed; as a consequence, the delta function of the Fermi's golden rule is 



replaced by broad spectral functions |15[]. We stress that this scenario, intimately related to 
the vector-potential or accelerated picture, has no counterpart in the scalar-potential one. 
Indeed, within the Wannier-Stark basis there is no carrier drift, and energy conservation is 



preserved. It is thus clear that such an effective semiclassical description of the ICFE is not 



gauge invariant [ |T6|j . 

Aim of this paper is to explain and remove this apparent contraddiction by providing a 
rigorous — i.e., gauge-invariant — formulation of Fermi's golden rule. Our analysis will show 
that the ICFE, as usually accounted for within semiclassical Monte Carlo simulations [IT5 



does not exist: it is simply an artifact due to the neglect of the time variation of the basis 
states which, in turn, leads to a ill-defined Markov limit in the carrier-phonon interaction 
process. This may account for the surprisingly good agreement between semiclassical and 
rigorous quantum-transport calculations |T0| . [T2| . 

In order to properly describe high-field transport in semiconductors, an electron-phonon 
system can be considered, whose Hamiltonian can be schematically written as 

H = H + H' = (H c + Hp) + H cp . (1) 

Here, the single-particle term H D includes the free-carrier and phonon Hamiltonians while 
the many-body contribution H' accounts for carrier-phonon coupling. More specifically, the 
free-carrier Hamiltonian 

H c = ^ ") + e ^(r, t; V ) + V\v) (2) 

describes the non- interacting carrier system within the periodic crystal potential V l (r) in 
the presence of a homogeneous external field F, the latter being included in a fully gauge- 
invariant form through the following vector and scalar potentials: 

A(r,i;r/) = -cqFt , ip(r; rj) = (77 - 1)F ■ r . (3) 

Here, the gauge freedom is expressed in terms of the transformation parameter r/. Indeed, 
it is easy to verify that for any value of rj — corresponding to a different Hamiltonian H c in 
Eq. (0) — we are describing exactly the same electric field F. In particular, for 77 = and 
77 = 1 one recovers the conventional scalar- and vector-potential formulations, respectively, 
which in turn correspond to the well-known Wannier-Stark and Bloch-oscillation pictures 



Let us finally introduce the carrier-phonon interaction Hamiltonian: 



H cp = E7 q b^ + b\e-^ , (4) 



q 



where the operators 6 q (6 q ) describe the creation (destruction) of phonons with wavevector 



q 



In order to provide a gauge-invariant formulation of carrier-phonon scattering, let us 
consider as basis states the (gauge-dependent) eigenstates of the free-carrier Hamiltonian: 



Here, the quantum numbers a — and therefore the corresponding eigenf unctions a (r) = 
(r\a) and energies e a — are functions of the transformation parameter rj and in general (i.e., 
for 7] 7^ 0) are also functions of time. Contrary to the conventional time-dependent pertur- 
bation theory, we thus propose a kinetic description based on a time-dependent quantum- 
mechanical representation. In particular, for rj = (scalar-potential gauge) we recover the 
well-known Wannier-Stark ladder |J: e a = e n = €q + nAe with Ae = eFd, d denoting the 
crystal periodicity along the field direction; In contrast, for rj = 1 we deal with the Houston 
or accelerated Bloch states ||: e a — e^), where k(£) = k + kt is the instantaneous carrier 
wavevector, k = eF/h being its field- induced time variation. 

Generally speaking, we shall denote with tP^ the unitary transformation connecting the 
eigenstates \a) in different gauges: 



Given such basis states {|a)}, most of the physical quantities we are interested in — e.g., 
carrier drift velocity and mean kinetic energy — are properly described by the single-particle 
density matrix [ |T7|j 



where (V a (a a ) denote creation (destruction) operators for a carrier in state a. This is 
defined as the average value of two creation and destruction operators: its diagonal elements 
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(5) 



a ( v )) =W' n \a{rj)) . 



(6) 




(7) 



fa = Paa correspond to the usual distribution functions of the semiclassical Boltzmann 
theory while the off-diagonal terms (a ^ a') describe the degree of quantum-mechanical 
phase coherence between states a and a'. It is easy to show that the density matrix (|7|) will 
gauge transform according to: 

Pa 1 a 2 = ^aia 3 Ps 3 a4^a4a2 > (8) 

where U-^a' — (ot\a') = ^a\U^ ,r, \a') are the matrix elements of U^ ,ri in the r\ representation. 
Here, the compact notation a = a(fj) has been introduced. 

For a time-dependent basis set {|ot)}, the Heisenberg equations of motion for the opera- 
tors a a are of the form ||: 

d d d 



dt a ° ~ dt ac 



* + dt aa 



(9) 



Compared to the standard equations of motion, the possible time variation of our basis 
states (p a gives rise to an additional term; the latter has been usually neglected, giving rise 
to the apparent discrepancies between scalar- and vector-potential gauges mentioned in the 
introductory part of the paper (see below). 

By combining Eqs. (0) and and considering the explicit form of the total Hamiltonian 
(HI), we get the following equation of motion for p: 



d d 



dt dt 



d 



Paa' — —i^aa'Paa' + ~7iPaa' +-J-Paa' (10) 



H CP dt 



with u aa i = (e Q — e a >) /%. The first, Liouville-like, term is due to the single-particle Hamil- 
tonian H D while the last one is again due to the possible time variation of the basis states 
a. The carrier-phonon contribution involves two-body as well as various phonon-assisted 
density matrices, e.g., s aa >^ = (a^6 q a a /) fL7f . These quantities describe many-particle cor- 
relations between carriers and phonons. Equation fllQ|) is thus the starting point of an infinite 
hierarchy involving higher-order density matrices. To obtain a solution — i.e., a closed set 



of equations — this hierarchy has to be truncated at some level. As clearly discussed in |T7 
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in order to properly describe carrier-phonon scattering, the time evolution of the phonon- 
assisted density matrix s aa / jq should be explicitely considered; its equation of motion has 
again the form of Eq. fllCf), i.e., 



(11) 



with £laa',q = ^^aa' ± ^q, ^q being the phonon dispersion. By treating the carrier-phonon 
contribution y^a'a Yia - ^ ne standard mean-field approximation |l7j, Eqs. (]ToD and ([□]) con- 
stitute a closed set of equations for the kinetic variables p and s; they form the basis of the 



so-called carrier-phonon quantum kinetics [17 



The semiclassical limit is finally obtained via an "adiabatic elimination" |17| of the 
phonon-assisted density matrices s. This consists in a formal integration of Eq. (|TT| ) on 
which a Markov limit is performed. More specifically, by neglecting the 0-term in Eq. (|i~T|), 



i.e., the contribution due to the time variation of the basis states a, the final result is: 



(12) 



with 



D(n m , q )= r dte-^o n ^^ t,)dt ' . 
Jo 



(13) 



By inserting the above formal solution for s into the carrier-phonon contribution of Eq. ([3~0|) 
we finally get a closed equation of motion for the single-particle density matrix p. In the 
linear regime, i.e., \p aa '\ "C 1, the carrier-phonon contribution to the dynamics is of the 
form: 



where 



d 
It 



Pc 



Hcp pp> 



(14) 



±,q a" 



(15a) 



and 



r„« W = 7T Y,Y1 S a'P9la",c l 9p'a",^' D (^t'c l ",n) ( 15b ) 
±,q a" 



are generalized in- and out-scattering rates [[17]. Here, the ± sign refers to phonon emission 



and absorption, respectively, J\f^ = iV q + \ ± \ denote the corresponding phonon occupation 
factors, and g a p^ are the matrix elements of the carrier-phonon Hamiltonian (|J). 

Equation ( |14D is the desired quantum-mechanical generalization of the well-known Boltz- 
mann transport equation ||14j| ; indeed, by neglecting all non-diagonal terms of the single- 
particle density matrix (p aa / = f a S aa i), the latter is easily recovered: 

jJ a \ Hcp = ^ ( T °0ff> ~ T pM ■ (16) 

Here, as usual, the scattering rates for in- and out-scattering processes coincide; they cor- 
respond to twice the diagonal parts (a(3 = a' f3') of the scattering operators P n and r out in 
Eq. 

r> = r> Q = 2tt]T \ 9aPA \ 2 ^ [v(n% f j\ . (17) 

The above semiclassical rates exhibit the well-known structure of the Fermi's golden rule; 
they describe the scattering probability for a phonon-induced transition between states a 
and f3. Their quantum-mechanical — or non-diagonal — generalization is then given by the 



scattering matrices (|15[) , which describe the effect on the time evolution of the density-matrix 
element p aa > due to the generic element ppp. 

The generalized carrier-phonon scattering rates in (|I5D — as well as their semiclassical 
counterparts in (|17]) — involve the D function in (0). For the case of a time- independent ba- 
sis set, i.e., i] = (Wannier-Stark states), the detuning frequency f2 is also time-independent 
and the real part of the function T> in (|13|) gives the well-known energy- conserving Dirac 
delta function of Fermi's golden rule; in contrast, for the case of a time-dependent basis, 
i.e., r] = 1 (accelerated Bloch states), the detuning is time-dependent, leading to a function 
T> with a Fresnel-like shape ||10|| . This is exactly the ICFE previously introduced ||: due 
to the field-induced variation of the carrier wavevector k, the energy difference between ini- 
tial and final states (ek(t) — Ck(t)±q) changes in time giving rise to broad resonances in the 



carrier-phonon scattering process. Such energy-nonconserving scenario has no counterpart 



in the Wannier-Stark picture [[0|. This clearly shows that the generalized scattering rates 
in (|T5| ) are not gauge invariant. 

Aim of this paper is to show that (i) the derivation recalled so far is only valid within the 
Wannier-Stark picture (77 = 0) and (ii) the ICFE previously described is simply an artifact 
due to the approximation scheme usually considered. Indeed, as anticipated, the crucial 
point is the neglect of the 0-terms, i.e., of the possible time variation of our basis states a. 
More specifically, a proper inclusion of these terms leads to a modified version of Eq. ([iTf): 

d 



with 

S aa ' being the matrix elements of the single-particle Hamiltonian (0) for rj — 0, i.e., written 
in the scalar-potential gauge. It follows that for a generic time-dependent basis, Eq. ( |T8"D has 
a non-diagonal structure, i.e., it does not allow a simple exponential solution. This implies 
that for 77 7^ the Markov limit is not straightforward. Indeed, the rigorous procedure 
is: (i) to perform a unitary transformation which diagonalizes the superoperator Q in (|19|). 
and (ii) to perform the exponential formal integration described above. Since S aa / are the 
matrix elements of H c for 77 = (scalar-potential gauge), the unitary transformation that 
diagonalizes Q is just U 0,v , i.e., the transformation connecting the generic gauge 77 to the 
scalar-potential basis (77 = 0). We stress that the new diagonal elements coincide with the 
eigenvalues of Q which, in turn, correspond to the time-independent detuning functions fl a a>,q 
in the Wannier-Stark gauge. This clearly shows that the Markov limit used to derive the 
generalized Boltzmann equation in ([H]) is only well-defined in the Wannier-Stark picture, for 
which the various 0-terms vanish and the detuning functions Q are time-independent. This 
does not violate the gauge-invariant nature of our formulation. Indeed, given the generalized 
Boltzmann equation fllil) written in the scalar-potential picture, the latter can be written 
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in any generic gauge rj by applying the unitary transformation U introduced in (|j). More 
specifically, let us introduce the single-particle density-matrix operator 

P = \ a )Paar(<x'\ , (20) 

act' 

which is by definition gauge invariant, i.e., //-independent. This can be easily checked by 
combining Eqs. © and @. This suggests to write the generalized Boltzmann equation (]T4|) 
in an operatorial form as: 

d 



dt P 



r> - r ou > + h.c. , (2i) 



where 



fWou t = ^ \a)\(3)T^(a'm (22) 

are in- and out-scattering superoperators. Our aim is to propose a gauge-invariant formula- 
tion of the problem. This requires the superoperators in fl2"2|) to be 77-independent as well. 



The analysis presented so far has shown that the scattering matrices f are only well 
defined in the Wannier-Stark picture (77 = 0). This allows us to extend their definitions to 
the generic 77 representation according to: 

where U-^a' are the matrix elements of the gauge transformation U 0,v in the Wannier-Stark 
picture (rj = 0). Here, implicit summation over repeated indices is assumed. 

Equation (^) is the gauge-invariant formulation of Fermi's golden rule we were looking 
for. Contrary to the conventional approach, in the case of a time-dependent basis, e.g., 
accelerated Bloch states Q, instead of using Eq. (|T5|) with an ad hoc energy- nonconserving 
T> function, the rigorous procedure is to compute the generalized scattering rates ( |I5| ) in 
the Wannier-Stark picture, and then to apply the gauge transformation U.^ according to 
Eq. 

In summary, we have proposed a gauge-invariant formulation of carrier-phonon interac- 
tion. It has been shown that ICFE — as usually described and accounted for — is just an 
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artefact due to the neglect of the time variation of the basis set. As anticipated, this may 
account for the surprisingly good agreement between semiclassical and rigorous quantum- 



transport calculations reported in [|1(J and |P2 |, as well as for the anomalous carrier heating 



typical of semiclassical ICFE models ||15|| . 

From our analysis we can conclude that the most severe approximation of the Boltzmann 
transport theory is not the Markov limit but the semiclassical approximation, i.e., the neglect 
of non-diagonal density- matrix elements. The latter, being intrinsically basis dependent, is 
not compatible with a gauge-invariant formulation of the problem. 
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